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Abstract 

The hierarchical reference theory (HRT) and the self-consistent Ornstein-Zernike approximation 
(SCOZA) are two liquid state theories that both furnish a largely satisfactory description of the 
critical region as well as the phase coexistence and equation of state in general. Furthermore, 
there are a number of similarities that suggest the possibility of a unification of both theories. 
Earlier in this respect we have studied consistency between the internal energy and free energy 
routes. As a next step toward this goal we here consider consistency with the compressibility route 
too, but we restrict explicit evaluations to a model whose exact solution is known showing that a 
unification works in that case. The model in question is the mean spherical model (MSM) which 
we here extend to a generalized MSM. For this case, we show that the correct solutions can be 
recovered from suitable boundary conditions through either of SCOZA or HRT alone as well as by 
the combined theory. Furthermore, the relation between the HRT-SCOZA equations and those of 
SCOZA and HRT becomes transparent. 
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I. INTRODUCTION 



Both the self-consistent Ornstein-Zernike approximation (SCOZA) [j], 0, 0] and the hier- 
archical reference theory (HRT) |4| have been found to give very accurate results for fluids 
in thermal equilibrium. In particular, the respective non-linear partial differential equations 
can be solved in the critical region, and their solution gives non-classical, and partly Ising- 
like, critical indices. These equations are derived by obtaining the equation of state in two 
independent ways and using thermodynamic consistency to fix a free parameter in the direct 
correlation function. 

Although both approaches appear similar in various ways, there are also marked dif- 
ferences. Both approaches make use of the compressibility route to thermodynamics, but 
SCOZA combines it with the internal energy route while HRT, inspired by momentum- 
space renormalization group theory, uses the Helmholtz free energy route. Thus, in short, 
the SCOZA adds effective strength to the attractive interaction by increasing inverse tem- 
perature (3=1/ kgT while HRT adds contributions to the interaction by including its Fourier 
components for shorter wavenumbers Q until the limit of interest Q — > is obtained. 

In a recent work we considered thermodynamic consistency between the internal energy 
and free energy routes to thermodynamics j5[. In the present work we want to extend 
this to consistency with the compressibility route, too. This requires the introduction of 
two free parameters instead of a single one in each of the original theories. In view of the 
high accuracy of HRT and SCOZA, one may expect this increased freedom to give even 
better results both for spin systems and for fluids. Due to the complexity of the combined 
problem, we here limit ourselves to a simpler situation that can be analysed explicitly and 
for which the exact solution can be established. This is the case for the mean spherical model 
(MSM). This model can be considered as the limit D — > oo for D-dimensional spins in d 
space dimensions where the transverse susceptibility is the relevant one for the fluctuation 
theorem or compressibility relation. In this connection we realize that the MSM can be 
extended in a straightforward way to a generalized MSM (GMSM) that yields the same 
HRT and SCOZA problems as the MSM; the only difference lies in the reference system 
boundary conditions. In the former the spin length is fixed to 1 while in the latter the 
spin length has some distribution of spin lengths. It should be pointed out that neither 
SCOZA nor HRT, nor our combined theory, are restricted to simple fluids and their lattice 
gas version although they are most often applied to these systems. The usual lattice gas is 
equivalent to the Ising model with spins s = ±1. So what we do here is to generalize and 
apply both of these theories to continuous spins on a lattice too. 

In Sec. |TT] we briefly consider the MSM, and in Sec. II I II we extend it to a GMSM whose 
solution is established. In Sec. HVI the SCOZA problem for the GMSM model is considered 
while in Sec. |V]the corresponding HRT problem is considered. With both of these approaches 
one free parameter can be determined. By appropriate choice of this parameter we get 
partial differential equations whose solutions are those of the GMSM with the reference 
system as boundary condition. In Sec. IVII an alternative method of solution based on a 
result from Ref. |5| is used. Then unification of SCOZA and HRT is considered in full 
generality in Sec. IVHI In Sec. IVIIII a pair correlation function of MSM form containing 
two free parameters is proposed and explicit HRT-SCOZA equations are established for the 
GMSM situation where transverse susceptibility replaces susceptibility. By analysis of these 
equations we show how the HRT-SCOZA equations work in this case and how the GMSM 
solution is recovered. 
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II. MEAN-SPHERICAL MODEL 



Consider D-dimensional spins on a lattice in d space dimensions and with cells of unit 
volume. It is well established that in the limit D — > oo this statistical mechanical system can 
be solved exactly 0, 0| • Its solution is the same as that of the spherical model [$] . Here we 
will consider its variant, the mean-spherical model (MSM) (9|. In the MSM the Ising spins 
are replaced by interacting spins s, whose length is Gaussian distributed such that (s?) = 1. 
This is nothing but a Gaussian model with an adjustable one-particle harmonic potential 
to keep (sf) = 1 fixed. More precisely one Laplace-transforms the spherical constraint by 
which a Gaussian model partition function is obtained. In the thermodynamic limit (i.e., 
for an infinite system) the inverse transform is determined by the maximum term of the 
integrand. Following the evaluations by H0ye and Stell the resulting Gibbs free energy g 



per spin becomes [10 



L = -(3g = s + + -lnvr--— ^- f In (2s - Pi>(k)) dk. (2.1) 

2(2s-/#(0)) 2 2(2vr) rf y 1 ^ " 1 ; 

Here s is the Laplace-transform variable, H is magnetic field, and ip{k) is the Fourier trans- 
form of the pair interaction normalized to ^(0) = 1. The maximum of expression (12. ip is 
obtained by taking dL/ds = 0. Then the equation of state follows easily by utilizing the 
condition for maximum. However, in the next section we will generalize the MSM so we will 
come back to the equation of state there. 



III. GENERALIZED MEAN-SPHERICAL MODEL. 

In the MSM the spherical constraint (sf) = 1 is fixed. For D-dimensional spins (D — > oo) 
this also corresponds to spins of fixed length. Now we can generalize this and let the 
D- dimensional spins have a distribution of lengths by which the spherical constraint is 
removed. At thermal equilibrium the average spin length squared will then change with 
both magnetization and temperature. This variation will depend upon the equation of state 
or the spin distribution specified for non-interacting spins. In the limit D — > oo this model 
will again be exactly solvable as the spin distribution becomes Gaussian, but the width of 
this distribution varies both with temperature and magnetization. 

With this in mind we first generalize the MSM to (sf) = n which replaces the first term 
of expression (12.11) by sn. Further, as n will not be fixed there must be a function F(n) that 
accounts for the distribution of n-values. In this way expression (12. ip can be generalized to 

L = sn + + -lnn--— i— f In (2s - (3${k)) dk - -F(n). (3.1) 

2(2s-/#(0)) 2 2{27i) d J 1 " 2 y ' y ' 

Again in the termodynamic limit the free energy is determined by maximum of expression 
(13. ip . but now with respect to both s and n. This gives the conditions 

dL ( f3H V If dk 



ds n \2s-p) (2-kY J 2s - (3${k) ° 

9 i = s- 1 -F>(n) = 0. (3.2) 
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With this the magnetization becomes 

m --?^-J*- (33) 

d({3H) ~2s-/3- {6 - 6) 

Now we put 

z = — and P(z) = — L= f i , (3.4) 

and Eqs. (13.21) and (13. 3p can be written 

P(n - m 2 ) = zP{z) 

| = fin) (3.5) 

with f{n) = 1/F'(n) and 

(3H = m(2s -f3)= m P(z) - (3m = m ( - - (3) . (3.6) 
n — m z V z } 



From this we obtain the transverse susceptibility x± as (f° r -D-dimensional spins) [10 

13 8(J3H±) (3H (3 



Xx dm^ m z 



(3 (3.7) 



where H± (— > 0) and m±_ (— > 0) are transverse magnetic field and transverse magnetization. 
The internal energy (from pair interactions) becomes 



fdL\ TT 11/" ip(k)dk 1 / 13 H 

U — — — — — mii = — 
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V«9/?; H 2(2n) d J 2s-(3^{k) 2\2s-P 

= --m 2 - —(P(z) - 1). (3.8) 

In accordance with this the spin correlation function for transverse correlations is 

f ± (Jfe) = . (3.9) 

(3(1 -zm) 

Further in accordance with the fluctuation theorem Tj_(0) = x±/P which is consistent with 
Eq. (13.71) . By integrating this using (13 Ah and (13.51) one finds the "core" condition 

r ± (0) = ^P(z) =n-m 2 . (3.10) 

P 

With n = 1 fixed one is back to the constraint of the usual MSM. 

For n not fixed the function f(n) must be defined by or related to the reference system 
at (3 — where P(z) = 1 as then z — > 0. So with Eqs. (I3.5p - (l3.7p and (13. lQj) we have for 
the reference system 

H = fi{m 2 ) = = n -m 2 = - = f(n). (3.11) 

Thus 

n = n{m 2 ) = /i(m 2 ) + m 2 . (3-12) 
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For general (3 Eq. (13.51) then means 



It = ^(m 2 ) = n - m 2 = ^P(z) = f(n)P(z). (3.13) 



Now we can put 
by which Eq. (13. 13[) becomes 



P(z) = l + 2J and Me =i = /( n ) 



n - m 2 = (1 + 2 J) fi e = f(n) + 2/x e J, 

n - m 2 e = f{n) = fi e (3.14) 

where 

m 2 e = m 2 + 2fi e J. (3.15) 
Comparing Eqs. (13. 14j) and (13. 1 5j) with Eqs. (13. lip and (13.121) one sees that this implies 

A^e = fie{fn 2 ) = n ( m l) — ml = M m e)- (3.16) 

IV. SCOZA FOR THE MODEL. 

In using the Self-consistent Ornstein-Zernike approximation (SCOZA) one assumes that 
the direct correlation function is the same as for long-range forces outside hard cores, i.e., 
—/3ip(k), and that there is an additional term that contains a free parameter determined 
from the consistency requirements of SCOZA. Traditionally this has been used to replace the 
temperature with an effective one, but here we will rely on the parameter z already intro- 
duced above. The SCOZA equation for D-dimensional spins (D — ► oo) connects transverse 
correlations and internal energy, i.e., with the substitution u = m 2 

« (A) » (§*) = L°JL = 2 ™ (4J) 

up V X± J op \ m J m dm ou 
Inserting from Eqs. (13. 7\ and (I3.8p the SCOZA equation becomes 



or (P'(z) = dP/dz) 



8(3 13 v 'du 
whose equations for the characteristics are 

d(3 (3 du dz 

~J = ~z 2 P'(z) = T' 

One solution of these equations is 

//._ = 

P 



^- Z 4p'(z)^ = z (4.2) 



(4.3) 



// e = - a = d. (4.4) 
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The other solution follows from 

du = -\P\z) dz = -dP'(z) dz 

as 

C x P{z) = C 2 -u. (4.5) 

Comparing with the solution of the generalized mean-spherical model (GMSM) above one 
finds that Eqs. (14. 4p and (14. 5 p are identical to the exact solution (13. 5p and (13.131) with 
constants of integration 

Ci = f(n) and C 2 = n. (4.6) 

Thus for a given reference system the solution of the SCOZA problem will reproduce the 
exact result (13.161) . 



V. HRT FOR THE MODEL. 



For the HRT we will also use expression (13.91) for the correlation function where z again 
is the free parameter. By adding interaction at wavevector k = Q for decreasing Q while 
keeping (3 constant, one obtains the equation \^ 



01 



— = AnCQ 2 In [1 - zip(Q)) with C = - 



1 1 



dQ 



2 (2tt) : 



(5.1) 



for space dimensionality d = 3. The / = —f3f = L — f3Hm where / is Helmholtz free energy 
per spin while L = —fig where g, which also appears on the left hand side of Eq. (13.11) . is 
the Gibbs free energy per spin. For the transverse susceptibility we then have (u 



m 







-P 



1 dl 



-2- 



.dl 



z m dm du 

With (15. ip inserted we get the HRT self-consistency equation 

_d_ 

dQ 

dQ \z) ~ ^dQ 



(5.2) 



z 



-^— = Att-2CQ 2 



-2— (— 

du \dQ 

j(g) dz 
1 - ztp(Q)du 



(5.3) 



Its equations for the characteristics are 



P 

Again the solution is 



dQ 



4tt • 2CQ 2 



1 - zif)(Q) 



du 



dz 



(5.4) 



CxPiz) 



--Co 



u 



where now in the integrand in integral (13.41) for P(z) = P(z,Q) the ip(k) is replaced by 
for < k < Q. Comparing one sees that this is nothing but Eqs. (14.41) and (14. 5p . 
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VI. ALTERNATIVE METHOD. 



By requiring consistency between Helmholtz free energy and internal energy when chang- 
ing both (3 and Q Reiner and H0ye obtained a more general solution beyond the one of the 
mean spherical approximation (MSA). They found [5( 



-n+1 



oo 

(Me = LI + AnRU ( 6 '^ 



n=l 



where K = J / Li e and P(z) = 1 + 2 J. The coefficients A n do not depend upon (3 and Q. Note 
that here the / does not contain the reference system and mean field terms. This expression 
will also hold in the present case when imposing consistency with the compressibility, but 
now the A n will depend upon the boundary condition at /3 = 0. We find by use of (16.11) 

f)T J^, 1 r)A 

W = /Hf„(m) - m = 8H 0im) - £ — l7 ^K»» (6.2) 

n=l 

where Hq is the reference system plus the mean field contributions. So for the transverse 
susceptibility (13.71) we get (/x e = z/(3) 



m ll e fJL 



ou / — ' n + 1 ou 



(6.3) 



This equation together with (16. ip determines Li e , i. e., the coefficients A n for given /i can be 
found by iteration by comparing equal powers of K. However, this problem can be trans- 
formed into the solution of a differential equation. So by choosing u and K as independent 
variables Eqs. (16.11) and (16.31) can be differentiated with respect to u and K respectively to 
obtain 

H 2 e dK du 1 } 

This partial differential equation has the solution 

fi e = d and 2C\K = (C 2 - d) - u. (6.5) 

With Ci(l + 2C X K) = C l P(z) this is again solution (jO]) and (Ol) . 



VII. COMBINED SCOZA AND HRT. 

In Ref . H consistency between free energy and internal energy was used to determine 
a single free parameter. This gave rise to first order partial differential equations whose 
properties were studied more closely. Now we will require thermodynamic consistency with 
the compressibility route, too, so that a second free parameter can be determined. 

Thus, to be general, consider a function \l/(/3, Q, m) which is determined via two free but 
unknown parameters 

z = z([3 } Q,m) } 

v = u(J3,Q,m). (7.1) 
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For the determination of z and u, the derivatives of \1/ are given by known functions of 

P, Q, m, z, and v as 

i&p = X = X(J3,Q,m,z,v), 

v&q = Y = Y{P,Q,m,z,v), (7.2) 
m" = Z = Z(P,Q,m,z,v). 

Here and below the subscripts mean partial derivatives with respect to (3 and Q etc. while 
the double prime means second derivative with respect to magnetization m. For the GMSM 
the latter is replaced by the first derivative with respect to u — m 2 . 
By differentiation with respect to j3 and Q we now get 

dftp = Xp dp + X Q dQ + X m dm + X z dz + X v du, 

d^ Q = Y/3 dp + Y Q dQ + Y m dm + Y z dz + Y u du, (7.3) 

d^>" = Zpdp + Zq dQ + Z m dm + Z z dz + Z v du, 

where subscripts indicate partial derivatives with respect to z and v. With three unknowns 
z and v the set of equations (17.2p represents a rather complex problem. We then note as 
m Ref.H that use of the identity d^> p/dQ = d^Q/dp will simplify this, and we first obtain 

X Q + X z z Q + X v v Q = Y{3 + Y z z/3 + Y v v p . (7.4) 

Further 

Zp + Z z Z/3 + Z v Vj3 = X", 

Z Q + Z z z Q + Z v v Q = Y" (7.5) 



i- (X" - z zZ ? - z p ) , 

±r (Y" - Z z z Q - Z Q ) . (7.6) 

These expressions can be used to substitute the vp and Vq in (I7.4p . and by some rearrange- 
ment the following equation is obtained 

\Z V X Z — X u Z z )zq — {Z U Y Z — Y v Z z )zp+ 
X U Y" - Y V X" + Z V (X Q - Yp) - Z Q X U + ZpY v = 0. (7.7) 

The previous one-parameter approximations can be recognized in Eq. (17.41) when v is 
considered constant. Then with z as the free parameter and — I — —fif where / is free 
energy per particle, one finds from Eq. ( 17. 6ft that up = is the SCOZA equation, uq = 
is the HRT equation, while the remaining (non-zero) terms of (17.41) give the consistency 
between the free energy and internal energy routes considered in Ref. 0. But in general 
these various consistencies give different z so (17.41) itself will not be solved by using one 
of these, except for the GMSM considered in this work. Since we know the exact solution 



or 



dm" 

~dQ 



"0 



s 



for the GMSM it should be possible to revover it directly from (17.41) and (17. 7p using a pair 
correlation function containing two free parameters. 

To obtain the resulting HRT-SCOZA equation the X" and Y" must be evaluated. In the 
usual case we then have 

X X m -\- X z z m ~\~ X u is m ^ 

X X mm -\- 2,X Tnz z m -\- 2X mv v m -\- X zz z m -\- 

2X zv Z m V m 3C vv V rn -\- X z Z mm -\- X v V mm (7.8) 

with similar expression for Y" with X replaced by Y . One notes that the v mm term will 
cancel when this is inserted in ( 17. 7p . Thus the resulting HRT-SCOZA equation becomes a 
second-order partial differential equation for z with coefficients that depend on v and its first- 
order derivatives. This can then be treated iteratively, by starting with some approximate 
u, solving for z, and updating v according to Eq. (17.61) . Note that v m = for m = 1/2 due to 
the symmetry of lattice gases if we identify v with fi e as we will do below. Its influence may 
therefore be only perturbing in the updating process and thus not crucial for the problem of 
performing a numerical solution. Anyway, here we will not try to pursue this question or try 
to analyse the properties of the general HRT-SCOZA equation any further. Instead we focus 
on the simplified situation with the GMSM to show how the HRT-SCOZA equation solves 
this problem. As mentioned earlier the susceptibility is then replaced by the transverse 
susceptibility. As in Sees. [TV] and [V] we then put u = m 2 to get 

X" ~ * ^llu = ^\XzZu + X v v u + X u ] 

Y" -+2% = 2[Y z z u + Y v v u + Y u }. (7.9) 

For (17. 7p this amounts to the substitution 

X V Y'' - Y V X" - 2[{X V Y Z ~ Y v X z )z u + X V Y U - Y V X U ], (7.10) 

where now the v u term cancels. Thus we are left with a first order partial differential equation 
for z with free variables (3, Q and u. But v, that is determined via ( 17. 6p . is still present in 
the coefficients of (17. 7p . 



VIII. TWO-PARAMETER PAIR CORRELATION FUNCTION. 

A simple way to introduce two parameters in the correlation function is to modify 
Eq. (ELS) into 

r x (k) = — . v ; - — = = — 8.i 

l-H(k)P1>(k) l-zip(k) 

for k > Q and k = (i.e., r_|_(A;) = for < k < Q). This means that the "self-energy" 
function is (for all k) 

t(k) = = . (8.2) 

This assumed form of the correlation function for continuum spins with two free parameters 
v and z can also be used for continuum fluids and their lattice gas version too. Thus various 
HRT-SCOZA expressions we derive in this section are also valid in the latter cases before 
we again specialize to the GMSA below Eq. (I8.14p . 
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An interesting feature of expression (18.11) is the adjustable amplitude v to which the 
internal energy is proportional. This may influence critical properties. For SCOZA there 
is a generalized kind of scaling [11]. The independence of v from z may change this. Note 
that here the v is not tied to a core condition which we here omit for simplicity. Such an 
omission may not be crucial for qualitative properties. Anyway, at least for SCOZA itself, 
the core condition is not crucial in this respect [12]. 

With r_|_(fc) given above we can now evaluate the quantities that enter the HRT-SCOZA 
equation. With \P = / = —j3f where / is Helmholtz free energy per particle we have [|| 

x = I = 7 '« + \" f (8 ' 3) 

dl ~ 
Y = ^ = 47rCQ 2 ln(l-S(g)MQ)) 

= 47rCQ 2 [ln (1 - z$(Q)) - In (1 - (z - v(3)^{Q))} (8.4) 

dl 1 
> — 

'du f ± (0) 



dl 1 1 — z , , s „. 
2— = = (u = m 2 ) (8.5) 



with, for given Q, 



J(z) = \p(z)-1) = C [ Z M dk. 

k>Q 



Here X = — U is a modification of expression (13. 8p for the internal energy U, the Y is 
a modification of expression (15.11) . while Z is the corresponding modification of expression 
(15. 2p . From this we obtain the partial derivatives 

Y p = -vL(Q,Az), 

Y z = -L{Q,z)+L(Q,Az), (8.6) 
Y v = -(3L(Q,Az), Y u = 



where Az = z — u/3 and 



L(Q , 2) = _iM£) =4 ^-i(« . (8.7) 
z 9Q l-ztp(Q) 



Further, with J'(z) = —dJ(z)/dz 



X Q = -uL(Q,z), 

X , = J-(m = + V W , (8.8) 



and finally 



dz \ z J z 2 z 

X v = -J(z), X u = -, 
z 2 



Z R = 0, Z = 0, 



Z, = -, Z u = 1-^. (8.9) 



v v 2 
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For the GMSM case where (17.101) is used for the X" and Y" terms we now can evaluate the 
coefficients of the HRT-SCOZA equation (17. 7p to obtain 



Az Q - Bzp + 2Cz u + D = 



where the coefficients are 



A = Z V X Z - X V Z Z = A 1 + A 2 , 
Ai = — —J(z), 



vz^ 



A 2 = — J'(z) 



vz 



(8.10) 



(8.11) 



B = Z V Y Z -Y„Z Z = B 1 + B 2 , 
B x = ^—Al(Q, z), B 2 = l^^L{Q, Az) 



With (17TT0|) we have 



C — X V Y Z — Y V X Z — C\ + Q 



2- 



Ci = -J{z) 

z 



-L(Q,z) + ^L(Q,Az) 



Finally, 



D 

D 1 
D 2 



C 2 = ^J'(z)L(Q,Az). 



Z U (X Q - Y p ) - Z Q X V + Z P Y V + 2[X V Y U - Y V X U ) 
Z U (X Q -Y (3 )-Y U = D 1 + D 2 , 

V 

1 - Az 



v 



L(Q, Az) = vB 2 . 



(8.12) 



(8.13) 



(8.14) 



Now in the GMSM the the solution to be expected yields v = z/(3. This suggests to 
replace z with a new variable fi = z//3, which simplifies the remaining analysis since then 
the D-terms will join 5-terms, and Eq. (I8.10p becomes 



with 



E x + E 2 + E 3 = 



1 -z, 



(8.15) 



E x = A 2 fi Q - Bxnp = (3 — —[J'{z)fj,Q + f3L(Q, z)^ ] 



E 2 = Ax/jlq + 2C x \i u = —J(z) 



^H Q + 2L{Q,z)ji u 



-B 2f M P + 2C 2l i u = — MQ,Az)[-(3 2 ^ + 2z 1 J'{z) i i u 



(8.16) 



Now one notes that E 3 = is the SCOZA equation (14.21) as P(z) = 1 + 2 J and z = 
Likewise E 2 = is the HRT equation (15.31) . These equations have both the common GMSM 
solution given by (14.41) and (14. 5p . Noting further that 

1 zJ'(z) L(Q,z) 



l-Z 



J(z) 



L(Q,Az)~ 
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it follows that the GMSM solution also solves E\ = 0, and by that it solves the HRT-SCOZA 
equation (18 . 1 5[) too. Here it can be noted that E\ = is nothing but consistency between 
the internal energy and free energy routes investigated in Ref. |5|- One should expect that 
Eq. (18.151) have other solutions too with a third constant of integration C 3 besides the two 
in Eq. (14. 5p . But in any case the GMSM solution is sufficient here as it can be adjusted into 
the reference system boundary conditions (13.111) or (14.61) . 

Finally we have to show that z = (3v fulfills Eq. ( 17. 6p too. With /j = v = zj (3 and 
Eqs. (ESD and ([Q]) - (^|) we get 

VP = -^-[2(X z z u + X v \i u + X u ) - Z z Zfi - Zp] = — [2/i 2 J'(z)/i M - zfip)} 
Z u 1 — z 

Hp = 2fi 2 J'(z)fi u (8.17) 

and 

= -^-[2{Y z z u + Y u /i u + Y u ) - Z z z Q - Z Q ] = -^—[-2L(Q, z)fi u - -/xq)] 
fi Q = -2fizL(Q,z)fi u . (8.18) 

And Eqs. ( 18T7D and (18J8D are nothing but the SCOZA and HRT equations E 3 = and 
E2 = as given by (I8.16p . Thus altogether we have shown in detail how the GMSM solves 
the unified HRT-SCOZA equations. We have reason to believe that this demonstration of a 
model that can be solved exactly will be useful for the possible solution of the HRT-SCOZA 
problem more generally where the second derivatives of Eq. (17.81) should be used. Also other 
assumptions for the correlation function different from the simple expression (18.11) may then 
be useful or may be needed. 



IX. CONCLUSIONS 

In the present work general equations for the unified HRT-SCOZA problem have been 
established using a simple form of the pair correlation function containing two free param- 
eters. To analyse the problem in more detail we have considered an exactly solvable model, 
the MSM and its extension the GMSM that we introduce. This generalization is also natural 
in so far as the GMSM is merely a more general solution of the same HRT-SCOZA equa- 
tions. The reference system boundary conditions determines the resulting solution. The 
SCOZA and HRT problems for the GMSM are first considered separately, and then they 
are combined. By analysis of the unified HRT-SCOZA it is shown how it can reproduce the 
known exact solution of the GMSM: Given correct boundary conditions and a suitable pa- 
rameterization of the correlation function, HRT-SCOZA successfully traces the evolution of 
the free parameters. We expect the analysis of how the HRT-SCOZA works for this special 
case to be useful for the more general situation of possibly solving the unified HRT-SCOZA 
problem. 
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